Hard-core bosons in one-dimensional interacting topological bands 
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We study the hard-core bosons in one-dimensional (ID) interacting topological bands at different 
filling factors using exact diagonalization. At the filling factor v = 1 and in the presence of on-site 
Hubbard interaction, we find no sign of the existence of the bosonic topological phase, which is in 
contrast to the fermionic case. Instead by studying the momentum distribution and the condensate 
fraction we find a superfluid (SF) to Mott-insulator transition driven by the Hubbard interaction. 
At the filling factor v = 1/3 and in the presence of longer-ranged interactions, we identify the 
bosonic fractional topological phase (FTP) whose ground-states are characterized by the three-fold 
degeneracy and quantized total Berry phase, which is very similar to the fermionic case. Finally 
we discuss the reason of the different behaviors of hard-core bosons at different filling factors by 
mapping them to spinless fermions. Our these results can be realized in cold-atom experiments. 



I. INTRODUCTION 

Recently topological insulators (TIs) attract intense 
theoretical and experimental studie»3J2J. Till now many 
materials are found to be TIs. The properties of non- 
interacting TIs have been well understood and some of 
their important properties are verified by experiments^^Sl 
Meanwhile the effects of interactions in TIs begin to be 
explored numerically and analyticallySHl21 p^^ ^jj^g mean- 
field level, the interaction can be decoupled to gener- 
ate spin-orbit coupling and topological Mott insulator is 
realizecP. Numerical simulations using different meth- 
ods obtain consistent resultpK^l. The interacting topo- 
logical invariant is developed using Green's function and 
simplified formula is proposed in terms of the Green's 
function at zero frequency or in the presence of inver- 
sion symmetrylSHlSl^ ^]-^g effects of interactions on the 
topolo gical c lassification of free fermion systems are also 
studiecpmi. 

By analogy with TIs, models that exhibit nearly flat- 
band with non-trivial topology are constructed in differ- 
ent systems, in which fractional Chern insulators (FCIs) 
may be realized in the absence of external magnetic 
ggjj^iMSll The phase is characterized by a multi-fold de- 
generate ground-states with quantized total Chern num- 
ber. By combining the two copies of FCIs formed by spin 
up and down electrons, fractional TIs with time-reversal 
symmetry can be constructed, which will be another new 
quantum state of matter. 

In real materials, the properties are usually exhibited 
by electrons which are fermions, so most of the above 
studies are for fermions. It is also interesting to ask 
whether there exist similar topological phases in bosonic 
systems. The studies in two dimensions have appeared 
and the properties of har d-cor e bosons in the topologi- 
cal bands are investigatecPEll. in this paper, based on 
our study on the ID interacting fermionic modeP^El]^ ^q 
study the behavior of hard-core bosons in ID topological 
bands using exact diagonalization. 

We study the cases of the filling factor v ~ 1 and 
V = 1/3. For the case of i' = 1, we consider on-site 
Hubbard interaction. By calculating the energies of the 



lowest states, the Berry phase and the fidelity metric of 
the ground-states, we find no sign of the existence of 
the bosonic topological phase. We further calculate the 
momentum distribution and the condensate fraction and 
find a SF to Mott-insulator transition driven by the Hub- 
bard interaction. For the case of j^ = 1/3, we consider 
nearest-neighboring (NN) and next-nearest-neighboring 
(NNN) interactions. We identify the bosonic FTP whose 
ground-states are characterized by the three-fold degen- 
eracy and quantized total Berry phase. The obtained 
phase diagram is very similar to that of the correspond- 
ing fermionic system except the different critical values. 
Finally we discuss the reason of the different behaviors 
of hard-core bosons at different filling factors. 



II. THE MODEL 

Our starting point is the ID interacting tight-binding 
model filled with hard-core bosons'^'^, 

H = ^(Af-|-2B)*la,*,-^B*Ja,*,+£ (1) 

i i,x 



where a^, (Jz are Pauli matrices, ^i = (ci, di)'^ with Ci{di) 
hard-core boson annihilating operator at the site r^ and 
n'i{nf) is the number operator of the orbit c{d). In the 
fermionic version of the non-interacting model (C/ = 0), 
depending on the values of the parameters A, B and M 
the system can be a trivial insulator or a non-trivial in- 
sulator at half-filling. Though like spinless fermion the 
occuping number of hard-core bosons is or 1 per or- 
bit on each site, the hard-core bosons obey commuta- 
tion relation [c^ , c W [di , dM ) = at sites i ^ j but anti- 
commutation {ci,cy}{{di,d\}) = 1 on sole site i, which 
makes the hard-core bosons exhibit different properties 
from the fermions. In the following calculations, we focus 
on the parameters' region where the Hamiltonian Eq.(|l]) 
at f7 = has non-trivial fermionic topological phase and 




FIG. 1: (Color online) (a) The energies of the ground- and 
first-excited states vs U. (b) The Berry phase and the fidelity 
metric vs U. The results in bosonic (black) and fermionic 
(red) systems are compared. The parameters are A — B — 1, 
M — —1 and the system size is L = 8. 
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FIG. 3: Finite-size scaling of the condensate fraction /c at 
f/ = 0, 10. The parameters are the same with those in Figfl] 
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FIG. 2: (Color online)The momentum distributions at (a) 
t/ = and (b) (7 = 10 on different sizes. The parameters are 
the same with those in Figll] 



study the properties of hard-core bosons in the interact- 
ing topological bands at different fillings. 



III. THE FILLING FACTOR v = 1 

We first study the case of the filling factor v = 1 (we 
denote the number of particles as Np and the filling factor 
is i^ = Np/L). To characterize the possible phases and 
phase transitions in the system, we calculate the energies 
En of the two lowest states, the Berry phase 7 and the 
fidelity metric g of the ground-state. The Berry phase is 
defined as 7 = § ^{i^e\-^\ipe) with 6 the twisted bound- 
ary phase and its value 7 mod 2tt gets a nonzero value 
TT for topological phase while zero for trivial phas^^^HMl 
The fidility metric g is defined as g{V, SV) — ^ ^ (sv)^ — 
with the fidelity F{V,SV) = |(*o(^)l*o(^ + SV) the 
overlap of the two ground-state wave functions at V and 
V + SV"^ . When the topological band is filled with hard- 
core bosons, as U is increased, it shows in Fig[T] that 
the ground-state remains gapped and 7 = 0, g = all 
the way, indicating that no obvious phase transition hap- 
pens. This is in contrast to the fermionic case, where 
the Hubbard interaction U drives a topological phase 
transition^. 

So the topological property doesn't persist when hard- 



core bosons replace the fermions. To identify the bosonic 
phase, we study the momentum distribution, which is 
defined by the formula'^^ 



L-l 



ndk) = jYl 



(4c,+4d,)e''^*-^'\ 



with the momentum k = {2tt/L)1{1 ~ 0, 1, ..., L — 1) and 
the average (■ ■ ■ ) over the ground-state wave function. 
As has been known for free hopping bosons, the ground- 
state is a SF, which is characterized by the peak at zero- 
momentum state and its height strongly depending on L. 
In our case oiU = (see FigM, the momentum distribu- 
tion shows peaks at /c = 7r/2, 3tt/2 and theirs heights in- 
crease with the size L. Thus the results show that the sys- 
tem is in a trivial SF phase at v = 1. When the Hubbard 
interaction is turned on, the system is expected to ex- 
perience a phase transition to the Mott-insulator phase. 
It is from the result at C/ = 10 [see FigMb)] where the 
momentum distribution ni^{k) tends to be uniform and 
its values nearly don't change with L. To further char- 
acterize the phases, we measure the condensate fraction 
/c = (Ac + Aci)/Nb {Nb is the total number of the hard- 
core bosons) with A^, (A^j) the largest eigenvalue of the 
one-particle density matrix p'^j = {c\cj){pf^ = {d\dj)'^^. 
In FigjS] it shows that ai U — Q fc scales to a nonzero 
value in the thermodynamic limit while at C/ = 10 it 
scales to zero. So for small U the system has a nonzero 
SF density, while for large U the system becomes a Mott- 
insulator. We emphasize that the study of the SF to 
Mott-insulator transition needs scaling for systems with 
larger sizes, which is beyond the present method. 



IV. THE FILLING FACTOR v = 1/3 

Next we study the case of the filling factor v = 1/3. 
We drop the Hubbard interaction, but add NN and NNN 
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FIG. 4: (Color online) The phase diagram in the (Vi,V2) 
plane at z/ = 1/3 for different sizes. The inset shows the 
finite-size scaling of the critical value V^c at different Vi . Here 
A = B = 1 and M — —2 when the Hamiltonian Eq.(fl| has 
non-trivial flatband in the fermionic case. 



interactions to Eq.(ll|, which writes, 

(iJ) iiiJ)) 

where Ui = n1 -\- nf is the total number of hard-core 
bosons on site r^ and Vi, V2 are the strength of the 
interactions. We have carried out the calculations at 
1/ = 1/3, and find the bosonic FTP where the ground- 
state is three-fold degenerate. We first glance at the 
phase diagram in the (Vi,V2) plane, which is shown is 
Fig|4] By turning on Vi, the ground-state is three-fold 
degenerate and the bosonic FTP emerges. The ground- 
state is separated from higher eigenstates by a finite gap, 
whose value increases with the strength of Vi . After turn- 
ing on V2, the value of the gap is decreased and vanishes 
at a critical value V2C, which marks the boundary in the 
phase diagram. Finite-size scaling shows that the bosonic 
FTP exists in the thermodynamic limit (see the inset of 
Fig§. 

We note that the present phase diagram is very similar 
to that of the correspon ding fermionic system except the 
smaller critical values V2C • In FiglSta), it is shown more 
clearly: at small V2 the energies En of the two lowest 
states are almost the same in the two systems, while at 
larger V2 they are different. The main difference of hard- 
core bosons and fermions is the exchanging relation. So 
when the number of the particles is fewer, the exchanging 
between the particles is less possible and the hard-core 
bosons are more like spinless fermions. Also at small V2, 
there is one particle within each isolated bond, thus their 
difference is further weakened. It is interesting that the 
situation is similar in two dimensions where hard-core 
bosons in topological band at integer filling don't ex hibit 
topological phase, while at fractional fillings they dcP^. 

In momentum space, the degenerate ground-states are 
in different momentum sectors and are equally spaced 



uj 





(a) 


^^^:^ 


-5 




Y 


-6 


/ — 


— Hard-core boson 

— Fermion 



r—^ (b) .— . f 
\ j L=9 \ / 




3 



FIG. 5: (Color online) (a) The energies of the ground- and 
the first-excited states vs V2 at Vi = 3. (b) The ground-state 
energy of each momentum sector at Vi = 3 and V2 = 1.6. (c) 
The total Berry phase vs V2 at Vi = 3. The parameters are 
the same with those in Fig 4]except M = —1.999 in (c) when 
the band slightly departs the exact flatness. Here the system 
size is L = 9. 



with the interval of A'p, as shown in Figjslb). We also cal- 
culate the total Berry phase of the ground-states, which 
is shown in Figlsjc). It shows that the total Berry phase 
gets nontrivial value tt for small V2 and begins to be ran- 
dom between (0, tt) from a critical value Vic, at which the 
bosonic FTP is broken (here the randomness is due to the 
fact that the multi-fold degeneracy of the ground-states 
is greater than three). The obtained critical value Vic is 
in good consistent with that from the energy spectra. 



V. MAPPING HARD-CORE BOSONS TO 
SPINLESS FERMIONS 



For ID systems bosons and fermions can be trans- 
formed into each other and there has been example 
showing that the topological features can be manifested 
by Bose to Fermi statistics transmutations in other 
ID systemaSll. So in the following we map hard-core 
bosons to spinless fermions using Matsubara-Matsuda 
and Jordan- Wigner transformationaSSBol^ ^g gain some 
insights of the different behaviors of hard-core bosons at 
different filling factors from the mapped fermionic model. 
Using Matsubara-Matsuda transformation, the Hamilto- 
nian Eq.(fl]) can be mapped to a spin-1/2 one with the 
identification c\{d\) = 5,+ (5+), c,(d,) = ^-(^rJ and 
"■ic,d = 5 + 'S'fc li- The two-component system can be re- 
garded as a two-leg ladder with one component on each 
site^^'^. For two-leg ladders Jordan- Wigner transforma- 
tion can be applied directly when all sites are arranged in 
a ID sequence (the zigzag path in Fig|6|. Then we divide 
the ladder into two sublattices and introduce two species 
of spinless fermions a^ and /3i. The spin operators on the 
two sublattices transform as: 

5+ = ^tg»7r5:^.^.(ataj-H^]/3j)gi,ratai 




FIG. 6: Zigzag path in the two-leg ladder. 

Using the above transformation, besides the terms in 
Eq.Q the following additional terms containing 4— and 
6— fermion operators appear: 

AH = 2BJ2 a\a,+in,p - 25 ^ /?| A+in,+i„ (2) 

i i 

+2iA'^a\Pi+i{n^p + ti^+iq - 2ni+ian^p) + H.c. 

i 

with riia = a\ai and Uip = Pjpi- At low filhngs when 
there is no doubly occupying and neighboring, these addi- 
tional terms vanish, so hard-core bosons show the same 
behaviors with fcrmions. While at high fillings, these 
terms show their cfl^cct and answer for the absence of 
the topological properties at i^ = 1 in hard-core boson 
systems. 

VI. CONCLUSIONS 

We have studied the hard-core bosons in ID interacting 
topological bands at different filling factors. For the case 



of i^ = 1, we consider on-site Hubbard interaction. By 
calculating the energies of the lowest states, the Berry 
phase and the fidelity metric of the ground-states, we 
find no sign of the existence of the bosonic topological 
phase, which is in contrast to the fermionic case. To 
identify the phase of the ground-state, we further study 
the momentum distribution and the condensate fraction 
and find a SF to Mott-insulator transition driven by the 
Hubbard interaction. For the case oiv — 1/3, we add NN 
and NNN interaction instead. We identify the bosonic 
FTP whose ground-states are characterized by the three- 
fold degeneracy and quantized total Berry phase. We also 
find that the obtained phase diagram is very similar to 
that of the corresponding fermionic system except the 
different critical values. Finally we discuss the reason of 
the different behaviors of hard-core bosons at different 
filling factors. Though the model we study is artificial, 
due to the rapid development of the field of cold-atoms^, 
it is hopeful that the model is engineered and the phases 
it exhibits are studied experimentally. 
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